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ABSTRACT 

We present Spitzer high spectral resolution IRS spectroscopy of three positions in the carbon-rich 
outflow of post-AGB star HD 44179, better known as the Red Rectangle. Surprisingly, the spectra 
show some strong unknown mid-infrared resonances, in the 13-20 fxm range. The shape and position 
of these resonances varies with position in the nebula, and are not correlated with the PAH features. 
We conclude these features are due to oxygen-rich minerals, located in a region which is believed to 
be predominantly carbon-rich. We provide possible explanations for the presence of oxygen-rich dust 
in the carbon-rich outflows. Simple Mg-Fe-oxides (Mgi_2,Fe K 0) are suggested as carriers of these 
unidentified features. 

Subject headings: stars: individual (HD 44179) — stars: AGB and post-AGB — circumstellar matter 
— dust, extinction — ISM: lines and bands 
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1. INTRODUCTION 



Ey er after its initial discovery paper ijCohen et al.l 
11975ft . HD 44179, lovingly nick-named the Red Rectan- 
gle, has proven to be a truly enigmatic astrophysical 
object, and hosts many unique phenomena. HD 44179 
is a post- Asymptotic Giant Branch (AGB) star with a 
stellar companion and surrounded by a circumbinary 
disk and a biconical nebula. In the optical, HD 44179 
shows the spect ral signature of the broad extended red 
emiss ion (ERE; iSchmidt et alJll980l iVan Winckel et aD 
2002), of which the carrier is still unknown. Superposed 
on this are sharper emission features which may be the 
Diffuse Interstellar B ands normally seen in absorption 
ijScarrott et al.ll992|L In the infrared the nebula exhibits 
remarkable evidence of a mixed chemistry: spectroscopy 
obtained with the Infrared Space Observatory (ISO), 
shows the presence of both oxygen-rich dust, in particu- 
lar sili cates, and carbon-r ich components in the form of 
PAHs l|Waters et alll998D . The authors suggest that the 
oxygen-rich dust predominantly resides in the disk, while 
narrow band imaging shows that the PAH emission is 
extended and fol lows the outflows i|Bregman et al]ll993l 
iSloan et alll993lh The beautiful opti cal images obtained 
with the Hubble Space Telescope (iCohen et alJ l2004|) 
show that the morphology of the nebula is far from homo- 
geneous. Regularly spaced brightness enhancements ex- 
ist, perhaps related to the gas or dust densi ty in the out- 
flow. The ISO spectroscopy presented bv iWaters et al.l 
( 1998) lacks the spatial resolution to show any correla- 
tion between the density enhancements seen in the Hub- 
ble images and the chemical composition. In this work, 
we report Spitzer Space Telescope observations of the 
Red Rectangle, to study spatial chemical variations in 
the outflow. 
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2. OBSERVATIONS AND DATA REDUCTION 

2.1. Observations 

We used the In fraRed Spectrograph fIRS: lHouck et alJ 
I2004D on Spitzer l|Werner et al.l2004a|) to obtain high res- 
olution spectroscopy of some positions in the outflow of 
the Red Rectangle. The observed regions in the outflow 
are about 30" away from the central source, to avoid 
saturation by the central binary and its direct environ- 
ment. Here we report the first results obtained with the 
Short High (SH) module of the IRS in two lines-of-sight 
in the nebula; a full description of the Spitzer data will 
be presented in a forthcoming paper (Markwick-Kemper 
et al. in prep.). The first slit position probes the re- 
gion just outside the biconical X, while the second one 
probes inside the biconical region (Fig. 2J. The coor- 
dinates of these slit positions are a = 6/i19to59.94s, S 
= -10°37'56.0" (J2000) and a = 6/il9m59.18s, S = 
-10°37'50.5" (J2000), respectively. At both these po- 
sitions we observed 2 cycles of 30 seconds. The SH mod- 
ule covers the 10.0-19.5 /jm wavelength range, with a 
spectral resolution of R w 600. 

2.2. Data reduction and calibration 

The spectra were reduced using the SMART package 
ijHigdon et alJl200"H) . We performed our data reduction 
on the products of pipeline version S9.5. We worked 
with the droopres images, which is the least processed 
product available to the observer. The last few pro- 
cessing steps of the pipeline were repeated off-line and 
adjusted to the needs of our observations: First, we di- 
vided the droopres images by the flatfields, to get the 
equivalent of the Basic Calibrated Dataproduct (BCD). 
We applied SI 0.0 flatfields, which yielded less noisy re- 
sults, to the S9.5 droopres data; while no S10.0 droopres 
data were distributed to the observers. Subsequently, full 
aperture extraction using SMART was performed on these 
images, using the flux conversion and other calibration 
files provided by the Spitzer Science Center (SSC). Bad 
pixels were removed by hand. Finally, where the orders 
were overlapping in wavelength coverage we cut off the 
band edges according to the order optimalization pro- 
vided by the SSC. In addition to the full aperture extrac- 
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Fig. 1. — The slit positions of the IRS SH observations on an 
image of the Red Rectangle from the red band of the Digitized Sky 
Survey 2. The foreground image is tak en in the WFPC2 F622W 
band on the Hubble Space Telescope ICohen et al.ll200l) . The 
two slit positions discussed in the text, are indicated with the la- 
bels 'inside' and 'outside' with respect to the biconical morphology. 
Positions of the individual pixels in the spatial direction within the 
slit positions are indicated with numbers. The dashed box in the 
center i ndicates the aperture used with the ISO SWS observations 
IWaters et al.|[l99l) . 



tion, our data was sufficiently complex that performing 
sub-aperture extraction was necessary. This mode of ex- 
traction is not offered for the high-resolution modules in 
the standard smart version, however, modified the full- 
aperture extraction routine to extract a specified fraction 
of the aperture and produce a spectrum of the flux for 
that sub-aperture. We approximated the size of the two 
high resolution slits as five pixel widths in the spatial 
direction by two pixel widths in the dispersion direction. 
We then extracted five sub-apertures of an entire column 
in the dispersion direction with a width of one pixels in 
the spatial direction. The sum of the five sub-apertures 
is equal to the original full aperture extraction. 

Currently the IRS data calibration is optimized for 
point sources, and little attention has been given so far 
to the extraction of spectra from extended sources. Sev- 
eral problems arise when extended source calibration is 
attempted. First, no extended source flux conversion ta- 
bles are available at this time. Using the flux conversion 
tables for point sources affects the absolute flux levels to 
some extent, although we expect the relative flux levels 
to be reasonably accurate. Second, it is not possible to 
use the measurements of the standards stars to calibrate 
the flux levels of extended sources, since it will cause 
fringing. A two-dimensional spectral response function is 
necessary to overcome this problem, but is not currently 
available. Finally, geometric effects should be taken into 
account, since the slit projection on the detectors cov- 
ers only 2x5 pixels, and the point spread function is - 
depending on wavelength - about 3 pixels large, compa- 
rable in size to both the i ndivid ual pixels as well as the 
entire slit. iSmith et alJ l)2004l) discuss the wavelength 
dependence of the point spread function, causing an ad- 
ditional slope in the spectra. For the high-resolution 




12 14 16 18 12 14 16 18 
A (/Lim) 



Fig. 2. — The full aperture spectra obtained from the two slit 
positions, outside and inside of the biconical region respectively. 
Indica ted are the presen ce of the out-of-plane bending modes of 
PAHs IHonv et al.ll200H) . 



modules, the tilt of the slit with respect to the detec- 
tor combined with the small number of pixels in the slit 
causes additional low-frequency fringing in the spectra. 
While the first geometric effect may be easily corrected, 
the second effect is fundamentally uncorrectable, since 
the spatial distribution of photons on a pixel is unknown. 

Fig. [21 shows the full aperture extracted spectra of the 
two lines-of-sight of interest, inside and outside the bicon- 
ical region. The observed variation in spectral features 
between the two positions is very large, and variations 
occur even on the scale of the pixel size. In Fig. [21 we 
show how the spectra change within both apertures, by 
extracting 3 pixel wide sub-aperture spectra while step- 
ping along the slit one pixel at a time. 

3. SOLID STATE FEATURES IN THE SPECTRUM 
3.1. PAHs? 

It is straightforward to identify the features shortward 
of 13 /im present in all lines-of-sight with the out-of-plane 
bend ing modes of the polycyclic hydrocarbons (PAHs; 
e.g. IHonv et aI1l2001() . but these features are dwarfed 
by the spectral features detected at wavelengths longer 
than 13 /im. These new resonances show some interesting 
and unique characteristics. First of all, the bands are 
very broad (several microns) and are very strong with 
respect to the dust continuum. Moreover, there is not 
a single set of well defined wavelengths at which these 
resonances occur, but instead the peak positions seem to 
vary spatially. 

Interestingly, inspection of the ISO SWS spectroscopy 
shows no evidence for the presence of the newly detected 
features, perhaps because in ISO's large beam many 
lines-of-sight are averaged, and individual contributions 
at various wavelengths gets averaged out to a broad, shal- 
low plateau, or alternatively, the carrier of these features 
may only be present in detectable amounts at consider- 
able distance from the disk, a region not probed by ISO's 
line-of-sight (Fig.GJ. 

Here we will explore possible carriers for the new fea- 
tures. Although the outflow shows strong PAH emis- 
sion at 11.2 fim even at this large distance from the 
central star, we can rule out PAHs as the main car- 
rier. Emission from PAHs in the 13- 20 p,m wavelength 
range is known and is also variable l)Honv et al.ll2fl0il 
iVan Kerckhoven et alJ l20001. PAH emission from 15-20 
/im varies from a very broad plateau that perhaps can be 
disentangled into a number of completely separate res- 
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Fig. 3. — Sub-aperture extractions compared with opacities of 
Mg-Fe-oxides. The bottom and middle parts of the plot shows 3- 
pixel wide sub-aperture extractions obtained in the lines-of-sight 
outside and inside the biconical region respectively. Three sub- 
aperture extractions are presented, centered on pixel 2, 3, and 4. 
The spectra obtained from the position inside the biconical region 
are offset with 0.9 Jy for clarity. The top part shows Mie opac- 
ities of Mg-Fc-oxidcs (Henning ct al. 1995) in arbitrary units. In 
increasing wavelength of peak position and increasing Fe-content, 
they are the opacities of Mgi-^Fe^O, with x = 0; 0.4; 0.5; 0.7; 0.8; 

0. 9; and 1. Besides the PAH features around 11 /im, two broad res- 
onances can be seen in the spectra from inside the biconical region; 
a feature from 13-16 fim and a feature peaking longward of 19 fim. 
While the feature peaking longwards of 19 (im can be fitted with a 
few different Mg-Fe-oxides, the carrier of the 13-16 fim feature re- 
mains unidentified. The spectra in the bottom part show a strong 
feature peaking at ~ 17.5 /im, which originates mostly from pixel 

1, as can be deduced from comparison of the sub-aperture extrac- 
tions. The opacity of Mgo.6Feo.4O is overplotted. 



onances in the 15-20 /im range l|V an Kerc khoven et alJ 
l207)0t IWerner et alJI2004rl ISmith et alJI2004D . which do 
not vary in peak position. Although these PAH bands fall 
in the correct wavelength ran ge, they are much w eaker 
than the 11.2 /im PAH band IjPeeters et all 1200^1 . and 
their FWHMs do not match. The low intensities and the 
fixed wavelengths at which the components of the PAH 
plateau are found, we rule out PAHs as a significant con- 
tributor to the resonances observed. 

3.2. Mg-Fe-oxides 

The extreme breadth of the resonances observed in the 
13-20 /im range is characteristic for solid state carriers, 
and the shift in the peak positions can be explained by 
variations in the temperature, composition, grain shape, 
grain size or lattice structure. The strong feature-to- 
continuum ratio dictates that the carrier has clean res- 
onances without a strong continuum component, while 
the breadth implies an amorphous structure. Moreover, 
the opacities in the UV and optical need to be high 
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Fig. 4. — Change of peak position with changing composition 
of Mgi^^Fe^O. The pluses indicate the the peak position as a 
function of composition, with an average error bar of 0.15 fim. 
The number density ratio on the vertical axi s is the Fe-content x of 
the laboratory samples from Henning ct al. Il99ff). The solid line 
represents a least squares fit through these data points. The peak 
position outside the biconical region is measured and corresponds 
to an iron content of 26%, whereas the inside the biconical region, 
the peak positions are found longward of 19 fim (indicated with 
the arrow), suggesting an Fe content of 73% or more. 



enough so that the grains heat to sufficiently high tem- 
peratures, while the band strength in the material needs 
to be strong enough so that the amount of dust required 
to explain the feature does not violate abundance con- 
straints. 

Although it is believed that the outflows of the Red 
Rectangle are dominated by carbon-rich components, the 
features from 13-20 /jm show this picture is incomplete. 
Simple Mg-Fe-oxides are a possible carrier for these new 
features. In Fig. [3] the spectra from the region outside 
the biconical outflow are compared to a few Mg-Fe-oxide 
minerals calculated f rom optical properties measured by 
IHenning et al.l l|1995|) . Indeed, peak strength and width 
are easily fitted, while the peak position suggests that 
a composition slightly more Mg-rich than Mgo.6Feo.4O 
is likely. The situation inside the biconical region seems 
quite different. There are two clear differences; a distinct 
band centered near 15 /im and a strong rise starting near 
17 /im. Interestingly, features of both bands, such as 
peak position and onset of the rise, shift as a function of 
position in the slit, and the shift seems correlated. This 
is consistent with a changing Fe-content in a relatively 
Fe-rich oxide. Using the correlation between composi- 
tion and peak position for spherical grains (Fig. 0} it is 
possible to constrain more precisely the compositions of 
the oxides. We find that the oxides inside the biconical 
region have an Fe-content Mg<o.27Fe>o.730, while the 
oxides outside the biconical region have a more Mg-rich 
composition of Mgo.74Feo.26O. 

Taking the Mg-Fe oxides as an example, we determine 
that oxide grains with a grain size of 0.1 /im heat up 
to about 75 K at 30" of the centr al star, using = 
43_2s, Tcff = 7750 K, d = 710 pc l|Men'shchikov et alJ 
2002) and assuming no extinction toward the central star. 
To account for the observed flux levels in the features a 
dust mass of M dust w 2.4 x 10 28 g = 1.2 x 1O~ 5 M in 
the form of Mg-Fe-oxides is required in the aperture of 
each of the slit positions, an amoun t which seems large. 
The distance of 710 pc derived by iMen'shchikov et alJ 
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(I2002f) might be a n overestimation (see e.g. Uura et all 
119971 BHobbs et aU feOCM). which would affect the derived 
dust masses considerably. For a distance of 500 pc for 
instance, the dust can be heated to 150 K at 30", and 
an oxide mass of only ~ 9 x 10 25 g, which is comparable 
to the mass of our moon, is needed in each aperture to 
account for the observed emission features. 

3.3. The 13-17 [im band 

The strong feature at 13-17 /xm remains unidenti- 
fied. However, because the peak position shifts with 
the change in onset of the feature longward of 17 /im 
we conclude that the two carriers are related, and that 
perhaps the 13-17 fim feature is caused by a mixture 
of Mg-Fe-oxides with another material. A good candi- 
date are spinels (MgA LjO^, which h ave features in this 
wavelength range ( Fab ian et all200H) . Their UV/optical 
opacity is too low to provide enough grain heating to ex- 
plain the observed flux levels on their own, but when they 
are in thermal contact with oxides, the spinels may get 
warm enough to explain the observed flux level, while 
the changing Fe-content still explains the shift in peak 
position. 

4. DISCUSSION 

The detection of new mid-infrared features in the bipo- 
lar outflow of HD 44179 is surprising. These features 
are carried by oxygen-rich minerals, possibly including 
simple Mg-Fe-oxides. Previously it was argued that the 
nature of the chemistry in these outflows is carbon-rich, 
while the dust in t he circumbinary disk is oxygen-rich 
ijWaters et al.l IT998J) . Independently obtained spatially 
resolved ground-based observations of the PAH features 
in the infrared support this view (e. g. iKerr et al.lll999t 
ISong et al.ll2003l IMivata et alj r2004). Our work, along 
with the re cent suggestion o f small PAHs in the circumbi- 
nary disk ijViih et alJl2005|) . indicate that the chemical 
distribution to the Red Rectangle may be more compli- 
cated than previously thought; both the disk and the 
outflow are not strictly oxygen-rich or carbon-rich, and 
apparently mixed chemistry environments occur in both 
components of the syste m. There are clear ly different 
episodes in mass ejection llCohen et al.ll200"H) and these 
may be connected to changes in the chemical composi- 
tion. 
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The origin of the oxygen-rich dust in the outflow of HD 
44179 remains speculative at best. It is possible that the 
oxides found in the outflow are remnants of an earlier, cir- 
cumstellar mass loss phase, which is currently being over- 
taken by the bipolar outflow. Simple oxides are found in 
the low dens ity winds of semi-r egular variables and early 
type Miras l)Posch et al .112002^ . likely progenitors of a 
system like the Red Rectangle. These stars are known 
to develop a outflow velocity which is significantly lower 
than what is seen in later st ages of the mass loss phase 
(e.g. iKnapp 8z Morrislll985|) . An alternate explanation 
is that the oxygen-rich dust found in the bipolar outflow 
originate from the circumbinary disk. Interaction be- 
tween the stellar wind and the disk may have eroded the 
disk and the wind may have dragged some of the disk ma- 
terial along in the collimated outflow directions. In this 
interaction, the oxygen-r ich dust grains fro m the disk 
(predominantly silicates; iWaters et al.1 Il998j) may have 
partially or completely evaporated, and recondensed high 
above the disk. When recondensation occurs, the den- 
sity will have dropped considerably, and the formation 
of significant amounts of anything beyond the simplest 
condensates may have been difficult, which could explain 
why the oxides are dominating the spectroscopy. 

Of course, whether Mg-Fe-oxides are indeed the carri- 
ers of these newly discovered oxygen-rich dust features 
remains to be seen, and further studies of the spectral 
appearance and spatial distribution of these features is 
required to further establish the nature of the dust. 



FMK is grateful for her time spent at UCLA as a 
Spitzer Fellow, working with Mike Jura. His detailed 
knowledge of the Red Rectangle and his insightful com- 
ments were of crucial importance to this project. We 
also thank Lou Allamandola for sharing his thoughts on 
the origin of the spectral features discussed in this pa- 
per. We thank Dan Watson and Ben Sargent for aiding 
in developing IRS data reduction tools, Jan Cami for 
careful reading of the manuscript and Hans Van Winckel 
for making available the HST WFPC2 image of the Red 
Rectangle. Support for this work was provided by NASA 
through the Spitzer Fellowship Program, under award 
011 808-001; and by the National Research Council. 



Men'shchikov, A. B., Schertl, D., Tuthill, P. G., Weigelt, G., & 

Yungelson, L. R. 2002, A&A, 393, 867 
Miyata, T., Kataza, H., Okamoto, Y. K, et al. 2004, A&A, 415, 

179 

Peeters, E., Mattioda, A. L., Hudgins, D. M., & Allamandola, L. J. 

2004, ApJ, 617, L65 
Posch, T., Kerschbaum, F., Mutschke, H., Dorschner, J., & Jager, 

C. 2002, A&A, 393, L7 
Scarrott, S. M., Watkin, S., Miles, J. R., & Sarre, P. J. 1992. 

MNRAS, 255, IIP 
Schmidt, G. D., Cohen, M., & Margon, B. 1980, ApJ, 239, L133 
Sloan, G. C, Grasdalen, G. L., & Levan, P. D. 1993, ApJ, 409, 412 
Smith, J. D. T., Dale, D. A., Armus, L., et al. 2004, ApJS, 154, 

199 

Song, I.-O., Kerr, T. H., McCombie, J., & Sarre, P. J. 2003, 

MNRAS, 346, LI 
Van Kerckhoven, C, Hony, S., Peeters, E., et al. 2000, A&A, 357, 

1013 

Van Winckel, H., Cohen, M., & Gull, T. R. 2002, A&A, 390, 147 
Vijh, U. P., Witt, A. N., & Gordon, K. D. 2005, ApJ, 619, 368 



New dust components in the Red Rectangle 



5 



Waters, L. B. F. M., Waelkens, C, van Winckel, H., et al. 1998, Werner, M. W., Uchida, K. I., Sellgren, K., ct al. 2004b, ApJS, 154, 

Nature, 391, 868 309 
Werner, M. W., Roellig, T. L., Low, F. J., ct al. 2004a, ApJS, 154, 

1 



